Current responses and voltage fluctuations 
in Josephson-junction systems 
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We consider arrays of Josephson junctions as well as single junctions in both the classical and 
quantum-mechanical regimes, and examine the generalized (frequency-dependent) resistance, which 
describes the dynamic responses of such Josephson-junction systems to external currents. It is shown 
that the generalized resistance and the power spectrum of voltage fluctuations are related via the 
fluctuation-dissipation theorem. Implications of the obtained relations are also discussed in various 
experimental situations. 



PACS Numbers: 74.50.+r, 74.25. Nf, 74.40. +k 



There has been much interest in the dynamics of 
Josephson junctions and Josephson-junction ar- 
rays P), e.g., current-voltage characteristics, dynamic 
resistivity, and voltage fluctuations. Among these, the 
voltage fluctuations provide direct information about the 
dynamic correlations in equilibrium whereas the 

resistivity probes the response to external currents 
The latter is also closely related to the relaxation func- 
tion, which describes the relaxation behavior towards the 
equilibrium state. These two probes are therefore com- 
plementary to each other, and one may expect, in view of 
the general idea of the fluctuation-dissipation (FD) theo- 
rem, that there exists a FD relation between them. Nev- 
ertheless most existing studies have been devoted cither 
to the resistivity or to the voltage fluctuations, and the 
relation between the two has hardly been investigated. 
Here we thus make use of the linear response theory to de- 
rive the generalized frequency-dependent resistance, and 
examine the relation between the generalized resistance 
and the power spectrum of the voltage fluctuations in 
Josephson-junction systems. 

There are three energy scales in a Josephson-junction 
system: the Josephson coupling energy Ej = hlj/2\e\, 
the self-charging energy Eq = e 2 /2C*o, and the junction- 
charging energy Ec = e 2 /2C, where Ij is the Josephson 
critical current and Co and C are the self-capacitance 
and the junction capacitance, respectively. In case that 
the Josephson energy dominates (E ,Ec <C Ej), the sys- 
tem at finite temperatures can be described by classical 
Langevin-type equations of motion for the phases of the 
order parameter on superconducting grains. For a single 
junction, in particular, the equation of motion reads 



+ 7^ + sin0 = / + C(t), 



(1) 



where <f> is the phase difference between the two grains, 
the time has been rescaled in units of the inverse of 
the Josephson plasma frequency uj p = h~ 1 \/8EcEj, I 
is the external direct-current bias in units of Ij, and 
the damping parameter 7 is related to the shunt resis- 



tance R via 7 = tuj p /2\e\RIj. The noise current £(t) (in 
units of Ij) is characterized by zero mean and correlation 
( COKO'j > = 2jT5(t-t r ), where T is the temperature in 
units of Ej. 

We introduce the probability distribution function 
P(4>, v, t) and write the Fokker-Planck (FP) equation ||, 
which corresponds to the Langevin equation (l|): 



d_ 
di 



P{cj>,v,t) = L FP {cl>,v)P{cj),v,t), 



(2) 



where v = <f> is the voltage (in units of hu) p /2\e\) across 
the junction, and the FP operator Lpp is defined to be 



L FP ((f>,v) 



d_ 



d_ 

dv 



jv + sin (/> — ! + jT 



(3) 



Suppose that the system is disturbed by a time- 
dependent external current SI(t) (in units of Ij), which 
gives the additional term —SI(t)(d/dv) in the FP opera- 
tor. The resulting change in the average value of v takes 
the form 



S(v(t)) = / dt' G(t-t')SI(t'), (4) 
J — 00 

where the linear response function G(t) is given by Q 
G(t) = -0(f) / deb dvv exp(fL FP ) j-Peqi^ v) 

JO J-oc UV 

(5) 

with the equilibrium probability distribution P eq ((j), v). 
Note that G(t) describes the voltage response of the sys- 
tem to the external current and can be expressed in terms 
of the correlation function 

G(t) = O^T- 1 ( «(t)[«(O)-a(0(O), w(0))] ) 

= 6(t)T^ 1 C v (t), (6) 



where 

In the frequency space, Eq. (g) takes the simple form 

2TRex(o;) = S v (u), (8) 

where the generalized resistance x( w ) and the voltage 
power spectrum S v (u>) are defined to be the Fourier trans- 
forms of G(t) and C v (t), respectively. Note also that it 
is the real part of x( w ) (rather than the imaginary part) 
which characterizes the dissipation across the junction. 
Equation (^) thus comprises the FD relation in the resis- 
tively and capacitively shunted junction (RCSJ) system, 
connecting the generalized resistance (i.e., dissipation) 
with the voltage correlation function describing equilib- 
rium fluctuations. 

The above relation is to be compared with the stan- 
dard FD theorem 0, which is applicable to a Hamil- 
tonian system. Namely, when one can explicitly write 
P oq = Z~ 1 e~* I ° !t ' T with a temperature-independent ef- 
fective Hamiltonian H c g, a(<p,v) in Eq. (Q) simply re- 
duces to a constant, giving the standard FD relation. In 
the RCSJ system, however, the periodicity in <fi (|] urges 
the system to have a non-zero probability current even 
in the stationary state. As a result, the stationary state 
cannot be described by an effective Hamiltonian, and in 
general a(<f>, v) becomes a dynamical variable depending 
on (j> and v. Similar features have been pointed out in 
the system without the ifi term ||, where the FD rela- 
tion between the correlation and response of <f> (rather 
than (j>) has been considered. Equation (||), in contrast, 
concerns the voltage v = (j>, which is a physical quantity. 
Here it is easy to show that the eigenfunction expansion 
of P oq |(| leads to the identity J dv (d/dv)P cq (<p,v) = 0, 
which implies that the (equilibrium) average of ot(cj>, v) is 
just the average voltage: ( a((f>,v) ) = (v) = v. It then 
follows that at long time scales the correlation function 
in Eq. (||) reduces to the standard voltage correlation 
function: C v (t) » C v {t) = ( [v(t)-v][v(Q)-v] ), which re- 
covers the standard FD theorem with the voltage power 
spectrum S v (uj) = dt e wt C v (t). 

On the other hand, in systems with the charging en- 
ergy dominant over the Josephson coupling energy, quan- 
tum fluctuation effects should be taken into account, es- 
pecially at very low temperatures. Such a macroscopic 
quantum system can be conveniently described by the 
quantum phase model: 



H = 4E c (n + q) 2 - Ej cos < 



(9) 



where the number n of excess Cooper pairs and the phase 
difference </> are quantum-mechanically conjugate vari- 
ables ( [n, <fi] = i), and q is the external gate charge in 
units of 2e. For simplicity, we assume that dissipation 
due to, e.g., quasiparticle tunneling is negligible in this 



low-temperature regime, and consider the system with- 
out the bias current, where q does not change with time. 
A disturbing current 8I(t) applied to the system leads 
to the perturbation Hamiltonian Hi = 8Ec(n + q) Sq(t), 
where 5q(t) = -K f dt' 61(f ) with K = y/Ej/8E c . 
The standard quantum theory of linear response then 
gives the induced voltage across the system in the form 



6{v(uj)) = Xq{u) Sq(w) = xM SI(u), 



(10) 



where x( w ) = (iK/u))x q (w) and Xg( w ) is the Fourier 
transform of the retarded Green's function G R (t — t') = 
i0{t - if) ( [v(t),v(t')] ) . The corresponding FD theorem 
reads 

(1 _ e -"/ KT )S v (u) =2Im Xq (u) = 2{u/K) Re X (w), 

(11) 

where the power spectrum S v (uj) is again given by the 
Fourier transform of the voltage correlation function 
C v (t), and lu should be understood as lim c ^ + (u + ie). It 
is pleasing that in the classical limit (oj/T — » 0), Eq. jll| ) 
reproduces the classical relation Eq. (g). Here without 
the bias current, we have v = and C~(t) — ( v(t)v(0) ). 
In the presence of the bias current, no increases with time 
and the unperturbed Hamiltonian Ho depends explicitly 
on time, which in general does not allow the standard 
derivation. Nevertheless when the bias current is suffi- 
ciently small, it may be incorporated into the disturbing 
current; this leads to the same FD relation as that shown 
inEq. (0). 

We now investigate the physical implications of the 
relations in Eqs. (||) and (|1J) to several cases. First, we 
consider the case that the external bias current is smaller 
than the Josephson critical current (I < 1). In view of 
Eq. (^), of particular interest in this case is the under- 
damped (7 -C 1) classical junctions at low temperatures. 
At zero temperature the phase of the (unperturbed) junc- 
tion stays at one of the local minima <f> = sin -1 / (mod 2tt) 
of the wash-board potential £/(</>) = — cos0 — I(j>. Small 
perturbations then induce the well-known plasma oscil- 
lation in the vicinity of the local minimum, with the os- 
cillation frequency u>p given by (1 — I 2 ) 1 / 4 at T = 7 = 0. 
At finite but sufficiently low temperatures (T <C 1), the 
noise current stirs up the small fluctuations around the 
local minimum, and the power spectrum of the voltage is 



known to become S v (u)) = 2 
@. The FD relation in Eq. (| 

RexM = 



Tw 2 [(w 2 - 
then gives 



>) 2 + 



7 2 u; 2 ] 



70; 



(w 2 



Ulp) 2 + ^ 2 UJ 2 '' 



(12) 



which reveals that at low temperatures the generalized 
resistance is temperature-independent. Note also the be- 
havior Rex(^) — > as u> — > 0, which is nothing but the 
manifestation of the superconducting channel. A system 
with strong quantum fluctuations displays another inter- 
esting phenomenon, the Bloch oscillation [rul, at suffi- 
ciently low temperature. In this case, Eq. (plj indicates 
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that like the voltage spectrum, the generalized resistance 
should also exhibit resonance peaks at the Bloch oscilla- 
tion frequencies proportional to the bias current. 

We next consider the Josephson oscillation in an over- 
damped (7 3> 1) Josephson junction with a large bias 
current (7 > 1) fujl , for which the Fokker-Planck equa- 
tion (g) reduces to || 

4 t p ^=^{ s ^- i+T ^) p ^- (i3) 

At zero temperature, the voltage v(t) across such a re- 
sistively shunted junction (RSJ) displays the Josephson 
oscillation with the frequency u>j = r )~ x \/l 2 — 1, which 
leads to the dc voltage v = u>j. As the temperature 
is increased from zero, the dc voltage also grows, ap- 
proaching the Ohmic characteristics. At the same time, 
thermal fluctuations introduce decaying behavior of the 
probability in addition to the oscillatory behavior. The 
eigenfunction expansion of the transition probability j(| 
at large time scales shows that the leading contribution 
comes from the lowest eigenvalue, giving 

P(4>, t\(f>o,Q) ~ e -7J ''' cos(ujjt) + higher harmonics, 

(14) 

where the line width 7j of the Josephson oscilla- 
tion are related to the dc voltage v(I, T) via 7,7 = 
^{dv/dlfT §[. While at T = 0, we obviously have 
7,7 = and us j = r y~ 1 ^/I 2 — 1, at finite temperatures, 
they may be estimated numerically to give 77 ss v = luj 
at T w 1 (see also Sec. 11.3 of Ref. §). Equation @ 
leads to the power spectrum 



jj + (w - UJj) 2 jj + {UJ + UJjf 

+ higher harmonics, 



(15) 



which shows correlation peaks at lu = ±luj; this is fol- 
lowed by resonances of the resistance via Eq. (||). It is 
also of interest to note the crossover behavior to the high- 
temperature (noise-dominant) regime at T sa 1, where 
the correlation peak in S v (u>) disappears. 

Heretofore we have concentrated on single-junction 
systems, but the generalization to an array system is 
straightforward. The classical equation of motion (Q) is 
easily generalized to any array of Josephson junctions, 
e.g., a two-dimensional (2D) square N x N array: 

j j j 

= 4 ( 16 ) 

with the noise current characterized by ( Qj (0) ) = 
2jT6(t)(SikSji — SuSjk), where Ay is the lattice Lapla- 
cian, dj = (Co/C)Sij + Ay is the dimensionless capaci- 
tance matrix, and the prime restricts the summation over 



the nearest neighbors of site i. Jj is the current fed into 
site i = (x,y), given by h = I(S xfi - 5 X , N ); this cor- 
responds to that along one edge of the array (x = 0), 
uniform current I is injected into each site, while along 
the opposite edge (x = N) the same current is extracted 
from each site. Similarly, in the quantum regime, the 
array is described by the Hamiltonian: 

v 

-Ej^coBfc-h) (17) 



<i]> 



with [rii, <f>j] = iSij, which is the obvious generalization 
of Eq. (g). With these, we can use the same procedure as 
that for a single junction, which leads to the conclusion 
that the FD relations in Ens. (||) and ( |TT| ) are applicable 
to junction arrays as well 12|. 

To examine the implications of the FD relations to ar- 
ray systems, we consider a 2D array at zero (direct) cur- 
rent bias, which is well known to display the Berezinskii- 
Kosterlitz-Thouless (BKT) transition Q at T = T B kt ■ 
According to the dynamic theory of the BKT transi- 
tion [|]|l4|] and ac electrical measurements |lq ], the con- 
tribution of the vortex bound pairs is screened out by 
the free vortices and the imaginary (inductive) part of 
the frequency-dependent complex impedance sharply de- 
creases to zero at frequency-dependent temperature T u 
(> Tbkt), which is accompanied by a peak of the real 
(resistive) part. The FD relation in Eq. (||) then suggests 
that the voltage power spectrum S v (u>) as a function of 
the temperature should also show a peak at T = T u . 
The voltage power spectrum, which, to our knowledge, 
has not been reported on a 2D array near the BKT tran- 
sition temperature, can be measured in equilibrium and 
complement another equilibrium measurement obtaining 
the flux noise spectrum ]l6[ ]. In the latter, an inter- 
esting relation between the flux noise spectrum and the 
frequency-dependent conductivity has been examined at 
T ~ Tbkt Wh - It would thus be an interesting topic 
in the future work to investigate in this region the con- 
nection between the voltage spectrum and the flux noise 
spectrum. 

For large bias currents (J > 1), the voltage spectrum 
S v (uj) has been studied numerically in 2D RSJ arrays [0, 
which indeed has revealed the correlation peaks as in 
Eq. (^5|). Further, it has been observed that the peaks 
in S v (oj) disappear at temperature T f» 1 (= Ej/ks), 
which happens to correspond to the BKT vortex unbind- 
ing transition temperature in the absence of the bias cur- 
rent. Here we point out that the apparent disappearance 
of the peak has little to do with the phase transition; 
rather Eq. ( |l5| ) shows that it is just a crossover. This is 
consistent with the result that there should be no phase 
transition in the 2D array driven by external currents 
larger than the junction critical current Eq]. 
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